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Abstract 

The minimization of transmission completion time for a given number of bits per user in an energy harvesting 
communication system, where energy harvesting instants are known in an offline manner is considered. An 
achievable rate region with structural properties satisfied by the 2-user AWGN Broadcast Channel capacity region 
is assumed. It is shown that even though all data are available at the beginning, a non-negative amount of energy 
from each energy harvest is deferred for later use such that the transmit power starts at its lowest value and rises as 
time progresses. The optimal scheduler ends the transmission to both users at the same time. Exploiting the special 
structure in the problem, the iterative offline algorithm, FlowRight, from earlier literature, is adapted and proved 
to solve this problem. The solution has polynomial complexity in the number of harvests used, and is observed to 
converge quickly on numerical examples. 

Index Terms 

Packet scheduling, energy harvesting, AWGN broadcast channel, flowright, energy-efficient scheduling. 

I. Introduction 

Since its formulation a decade ago [1], the problem of energy-efficient packet transmission scheduling 
has drawn considerable interest from the research community [0, flU. The basic offline problem 
is to assign transmission durations (equivalently, code rates) to a set of packets whose arrival times are 
known beforehand, so that they are all transmitted within a given time window with minimum total energy. 
Recently, the problem has been recast with a formulation where the goal is to minimize the time by which 
all packets are transmitted, given that energy is harvested at certain known instants flU. 

In this paper, we extend the formulation in [5] to a multiuser scenario with one sender and multiple 
receivers. In particular, we assume an AWGN Broadcast Channel where the sender gets replenished with 
arbitrary amounts of energy at arbitrary points in time. The harvested energy becomes instantly available 
for use, and the transmission power can be changed at any time by the sender. The choices of power 
level and the rates to individual receivers across time is called a schedule. The sender needs to transmit a 
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certain number of bits to each receiver. We consider the case that these bits are available at the beginning 
of transmission. The goal is to find a schedule that achieves the unique minimum time T opt , by which 
the data of all users can been transmitted using the given sequence of harvests. Throughout the paper, we 
focus on the offline problem, where the energy harvesting times as well as packet arrival times are known 
in advance. The online version of the problem in which the times of energy harvests are not known a 
priori and decisions need to be made in real-time as the harvests occur, is interesting yet analytically less 
tractable and left outside the scope of this paper. 

It is well known that both with optimal and practical coding schemes, the energy per bit increases with 
the transmission rate, in other words, transmitting fast is inefficient in terms of energy JS). This is the root 
of the sender's dilemma: it will pay off for the sender to slow down, yet it needs to minimize the overall 
transmission duration. Interestingly, it turns out that even if all packets were available in the beginning, 
the optimal schedule starts slowly, deferring some of the harvested energy for future use. More precisely, 
we will show that in the optimal schedule the transmission power is non-decreasing in time, similarly to 
the point-to-point schedule [0. 

In the point-to-point problem, determining power levels determines the schedule, as transmission rate 
is a function of average power. In the broadcast problem, however, there is no one-to-one correspondence 
between the transmission power and the rate point. For example, with optimal coding, there is a continuum 
of rates on the boundary of the capacity region corresponding to a certain average power constraint. Hence, 
the rates and the power have to be determined together. We observe that in the optimal schedule, the average 
rates used by the users are proportional to their numbers of bits, i.e. the schedule always continuously 
transmits to all users at the same time and finishes transmission to all users at the same time. Having 
made this observation, we can exploit the mathematical similarities between this problem and the problem 
in (6), and show that the solution is found by the algorithm FlowRight, defined in (6) and adapted here 
to work with different parameters. 

In the next section, we make observations about the two-user AWGN broadcast channel capacity region. 
The statement of the problem as a cost minimization problem, as well as its solution will use certain 
structural properties of the AWGN capacity region, such as the monotonicity and convexity of the average 
power with respect to the rate pair. Of course, this specific rate region can only be approached under 
optimal coding as blocklengths and the number of information bits go to infinity. For example, in the 
single user AWGN channel the numerical value of the minimum energy per bit corresponding to a given 



reliability monotonically decreases with the number of information bits UTOl . However, the basic structural 
properties of the rate region will be satisfied by the achievable rate regions of many suboptimal practical 
coding schemes as well as finite blocklength optimal coding schemes. 



We define the problem in Section III In Section IV we explore the properties of the optimal solution. 



This is followed by the description of the modified FlowRight algorithm, and the proof of its convergence 
and optimality of the resulting schedule. The complexity of the iterative algorithm is analyzed in Section 



|V} The implementation of this algorithm is discussed, followed by a numerical example in Section |VT 
Section IVTll summarizes our conclusions and outlines further directions. 

II. Broadcast Channel 

Consider a discrete-time AWGN broadcast channel with one sender and two receivers. The signal 
received by the i th user at time k is given by 

Y t [k] = ^X[k] + Z l [k], (1) 

where X[k] is the transmitted signal with average power constraint P, y^'s are the channel gains and the 
Zj[fc]'s are i.i.d. zero-mean Gaussian noise with variance a 2 . The capacity region of the channel assuming 
s\ and s 2 are constants and s\ > s 2 > , is the set of rate pairs (ri,r 2 ) such that 

ri < ilog 2 (l + ^) (2) 
1. (. (\- a)s 2 P \ 



r 2 < -log 2 l+ v ' \ ) (3) 
2 \ as 2 P -\- cr J 

for some < a < 1. Hence, the l s< user is the stronger user. 

It is straightforward to show that, given s 1 and s 2 , for any Pi > P 2 , the capacity region corresponding 
to an average power constraint Pi dominates the one corresponding to P 2 . Therefore, given a rate pair 
(r 1? r 2 ), there is a unique P = g(ri, r 2 ) (see fl6|) such that (r 1; r 2 ) lies on the boundary of the rate region 
with power constraint equal to P. After replacing the inequalities in (|2]) and (|3]) by equalities, the function 
g(ri,r 2 ) is written as follows. ([2]) can be written as as as\P/a 2 = 2 2ri — 1. Hence, a = (cr 2 /Psi)(2 2ri — 1). 
After substituting into ([2]) and rearranging the terms, we obtain 

\ S 2 Si J 

The function g(ri, r 2 ) is twice continuously differentiable and strictly convex in r x and r 2 . Throughout 



the paper, it will be useful to express the t\ and r 2 as a function of each other and the minimum power 
P. By algebraic manipulation of ([2]) and ([3]), we obtain the following: 

n-W 8 ) - -log 2 ( S2 a^ ^ (5) 

r 2 = h 2 (P, ri ) = ^log 2 ( f(2 |^_ + 1 1 ) + 1 ). (6) 

The properties satisfied by these rate functions for the AWGN BC capacity region with si > s 2 
summarized in the following will be used in the rest of the paper. 
Properties of the rate region: 

1) Nonnegativity: hx(P,r) > 0,h 2 (P,r) > 0. 

2) Monotonicity: hx(P, r), h 2 (P, r) are both monotone decreasing in r, and monotone increasing in P. 

3) Concavity: hi(P,r) and h 2 (P,r) are concave in P and r. 

4) The rate of the user with the weaker channel satisfies the following: ^dp^ = 0, = 0. 

Proposition 1: The functions hi and h 2 , defined in ([5]),(|6]) on x satisfy the Properties of the 
rate region given by (l)-(4). 
Proof. See Appendix A. 



III. PROBLEM DEFINITION 

Consider the broadcast link as described in the previous section, with a sender who needs to transmit 
P>i < oo and B 2 < oo bits with a certain degree of reliability to users 1 and 2, respectively Q Also 
assume that at time ti=0, sender has E\ > units of energy available and at times t 2 , ...,tk+i, energies 
are harvested with amounts E 2 , ...,E K+1 , respectively, as depicted in Fig [T] Inter- arrival times of energy 
harvests are named as epochs, and marked with = 1, k. 




Fig. 1. System model with all of the bits of the two users are available at the beginning. Energies arrive at times t, where i £ {1, 2, k + 1} 

'Throughout the paper, two receivers will be considered for ease of exposition. However, the results can be generalized to more than two 
receivers. 



It will be assumed that the sender has the ability to change its rate pair at any time, according to the 
available energy and remaining number of bits. Such ideal adaptation, which has been used in previous 
literature (e.g., [5], and references therein), may be approximated by using adaptive coding and modulation 
in a practical system. 

Starting at t — 0, let {(rn, r 2 i), (ri 2 , r 2 2), (Vin, r 2n ), ■ ■ ■}, be the successive pairs of rates used by 
the sender, and {h,l 2 , ■■■Jn, ■ ■ •} be the respective durations for which these pairs are used. Here, > 
is user i's rate in the j th rate pair. By definition, at least one user's rate changes from one rate pair to 
the next one. We will refer to the sequence of rate pairs and durations as a schedule. The problem of 
interest is to find an optimal offline schedule, that is, a schedule that minimizes the overall transmission 
completion time of the B\ and B 2 bits to their respective destinations, with complete knowledge of future 
energy harvesting instants and the amounts to be harvested. 

It will also be assumed that the problem is feasible; that is, sufficient energy will be harvested to 
transmit the given Bi < 00 and B 2 < 00 bits in arbitrarily large but finite total time, T. Note that for 
any given E, there is a small enough rate (equivalently, long enough transmission duration) such that 
B\ and B 2 bits can be transmitted with energy E, provided that the minimum energy per bit required 
for communication on the broadcast channel for the given finite amount of bits is satisfied |0. In the 
point-to-point case with infinite blocklengths, the well known limit for energy per bit is —1.59 dB. For 
sending finite amounts of data, the minimum energy per bit is higher even at nonvanishing values of error 
probability. However, the upper and lower bounds in iflOl on energy per bit come very close to the ideal 
limit at B = 10 3 bits, and even at smaller numbers of bits. 

In order to define the two-user broadcast channel offline scheduling problem as an optimization problem, 
we will use the set of observations stated in Lemmas [TJj5j Lemma [5] will establish that in an optimal 
schedule the transmission to both users ends at the same time. Lemma [2] will establish that in an optimal 
schedule the rates and power level do not change between energy harvests that are used. As a consequence 
of these two results, the problem reduces to Problem [TJ 

We start by proving a more general result than Lemma [2] which will be used in the proof of Lemma [2] 



as well as Theorem 1 in Section III Specifically, we take a finite time window which is divided into two 
slots such that different power levels are used in each. We show that by using a more even distribution 
of power (reducing the difference of the power levels) as much as energy causality permits, at least the 
same amount of data can be transmitted in the same amount of time using the same amount of energy. 
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In the special case when this time window is within (or all of) one epoch, all the energy that is used is 
available in the beginning hence the powers can be completely equalized. 

Lemma 1: Suppose that within a time window (ti, t 2 ), the sender changes its transmit power at point 
r* such that T\ < r* < t 2 . Keeping the total consumed energy in (ti,t 2 ) constant, the sender can send 
at least the same number of bits to the users within the same duration by bringing power levels closer to 
each other, if feasible (i.e., unless such a change requires energy to be used before its harvested.) 

Proof. Let the total duration be t = r 2 — T\, and the lengths of the two slots (3t and (1 — (3)t, with 

power levels in the two slots Pi and P 2 , as illustrated in Figj2j Denote the rate pairs in the l s< and 2 nd 

slots as (r n , r 2 i) and (r 12 , r 22 ), respectively. 

pt (i-p)t 
-< >- -< ■>- 

1 p h p 1 

X, P i X* p 2 x 2 

( r n> r 2i) (r 12 , r 22) 

Fig. 2. Illustration of the transmission scheme used in Lemma [T] 

First, consider the case where the power level used in the first slot is smaller: Pi < P 2 . When P\(3t is 
equal to the total energy available for use in r*), transferring energy from the second slot to the first 
is not feasible, and we stop. However, if it is possible to transfer some positive amount of energy AE 
from the second slot to the first, we shall show that we can only improve the allocation. 

Let us denote the average rates for the stronger and weaker users as f\ = (3rn + (1 — f3)r\ 2 and 
^2 — P r 2i + (1 — ft) r 22, respectively. We will show that keeping the total consumed energy and fi 
constant, the sender can achieve an average rate r 2 for the second user such that r 2 > r 2 by changing Pi 
to P[ and P 2 to P 2 satisfying 

Pi < P'x < P' 2 < P2- (7) 

While keeping the total energy constant, a certain amount of energy should be transferred from the 2 nd 
slot to the 1 st one in order to satisfy (|7]). In this case, we have the following 

P[ = P x + (1 - (3)AP ,P' 2 = P 2 - (3AP (8) 

Average rate belonging to the weaker user over the whole duration t is given by 

r = = / i2 (p;,ri)/3 + / i2 (P 2 ,ri)(l-/3) 

> h 2 (Pi,r n )P + h 2 (P 2 ,ry i )(l-P)=f2 



follows from the fact that 

h 2 (Pi n)/3 + h 2 (P 2 , fi)(l -13)- h 2 (P u r n )P ~ h 2 (P 2 , r 12 )(l - 0) > (9) 



for all (3 = {0, 1} with equality achieved at f3 = 0, 1. This is a consequence of the properties listed in 
Section [n} and proved in Section |IX-B in the Appendix. 



In the remaining case, P\ > P 2 , a similar argument holds where P 2 < P' 2 < P[ < P±. In this case it 
is always possible to strictly improve the allocation by transferring a positive amount of energy from the 
first slot to the second, as energy can always be deferred for future use. ■ 

Corollary 1: In a schedule that ends at T opt , power does not change within epochs in [0, T opt ]. 

Proof. The claim is that the power does not change within epochs, of course with the exception of the 
last epoch. (In the last epoch that is used, the transmission ends and the power is reduced to zero at some 
point within the epoch.) By definition, no new energy or data is added during an epoch, so it is intuitive 
that the decision on power allocation does not change at a point during an epoch. To reach contradiction, 
suppose that the sender changes its power allocation during an epoch. From Lemma [Tj the power levels 
can be allocated closer to each other so that at least the same number of bits can be transmitted to the 
users. Since this case is not limited by causality, this procedure can be continued until the power levels 
within the epoch are equalized, strictly improving the schedule, contradicting the optimality of the original 
schedule. ■ 

Lemma 2: In a schedule that ends at T opt , the rate pair does not change within epochs in [0, T opt ]. 

Proof. From Corollary [T] power level stays constant during epochs in an optimal schedule. Now, suppose 
the sender changes its rate pair at some point during an epoch, while the power is constant at P. Let the 
lengths of the two slots as (3t and (1 — (3)t and the rate pairs in the I s * and 2 nd slots as (rn, r 2 i) and 
( r i2> r 22)- Due to the concavity of h 2 (P,r) in r, setting r x to the average rate only improves r 2 

h 2 (P,/3r n + (1 - (3)r 12 ) > 0h 2 (P,r n ) + (1 - P)MP,r 12 ). 

Hence, by equating the rate pair, at least the same number of bits can be transmitted at the same time. ■ 
The next result is an observation of the structure of the basic solution when there is only one energy 

harvest (the one at t = 0). 

Lemma 3: Suppose E { = Vi > 1 in the system model in FigjT] To minimize the overall transmission 

duration, the sender finishes transmission to both users at the same time. 



Proof. To reach contradiction, suppose that in an optimal solution, the sender finishes transmission 
to one of the users before the other. This means that the rate pair changes at some point (when the 
transmission of one of the users ends before the other), although no new energy has been harvested. By 
Lemma [2j averaging the power levels and rates and using one rate pair continuously would enable us to 
send at least the same number of bits during the same time. This contradicts the optimality of the original 
solution. ■ 

Lemma [3] tells us that the ratio of the rates ri/r 2 is equal to the ratio of the bits Bi/B 2 . Then, for the 
AWGN case from (|2]) and (|3]) the ratio of powers, a, can be found by setting: 

ri = 2 l0g2(1 + ^ T2 = 2 l0g2(1 + ^T^ } 
where p 1 = Pa(Bi, B 2 ), p 2 = -P(l — cx(B^ B 2 )). Using r 1 /r 2 = Bi/B 2 , one can obtain 



1 + 



as\P 



a 1 



B 2 



1 — a 



Bi 



ols 2 P + a 2 J 

Solving for a and substituting into (|2]) and ([3]) yields a rate pair (r\, r 2 ), for a given value of P. 

Before discussing how to find the right value of P, it will be illustrative to present an alternative proof 
for Lemma [3] Suppose the sender has a power level P to use. The question is to obtain the minimum 
termination time for all the bits, T min , given by the following: 

Bi B 2 , 



T ■ 

-L mir 



nuii(max( — , — '-)) = f(ri,r 2 ) 
n r 2 

n •> r 2 B 2 
B2 if ri > Bi 




decreases 



Fig. 3. Illustration of Lemma[3] The rate pair (n, r 2 ) that minimizes overall transmission time for the sender to transmit Bi and B2 bits 
to each user is at point A.(cf. properties (l)-(4), the shaded rate region is convex.) 

The contours of constant f(r 1 ,r 2 ) are shown in Figure [3] The value of / gets smaller as we move 



outward from the rate region. The last contour that still touches the region, touches it at the point labelled 
A. Therefore, T min is obtained by the rate pair at point A on the boundary of the rate region, which 
satisfies r\jr 2 = B\jB 2 . 

Given E, B\, and B 2 , determining T min entails solving a nonlinear equation which can be done iteratively 
(see Algorithm [T]) using the bisection method. The total energy used to transmit B\ and B 2 bits, given by 
f(T) = T.g(r 1 ,r 2 ) = T.g(Bi/T, B 2 /T) is convex, monotonically decreasing in T. Combining this with 
our initial assumption about E being large enough to satisfy the minimum energy per bit requirement 
(E > liniT-K>o T.g(Bi/T, B 2 /T)), there is always a unique smallest value T for which T.g(ri,r 2 ) is just 
below E. In the bisection method, the objective is to find the single root of the equation f(T) — E = 0. 
Assume that we have an upper bound T uppcr > T min . Starting with the initial domain interval [T mm = 
, T max = T up per] j at each iteration the domain is bisected and the subinterval in which the root T min 
lies is selected as follows. If f(T) - E < 0, then we set T max = T, otherwise we set T min = T and 
T is selected to be the mid-point of the updated interval for the next iteration. Algorithm converges to 
the unique solution as the domain is continuous and can be terminated within a certain arbitrarily small 
tolerance e > in a practical implementation. This iterative method has been used in generating the 
numerical examples given later in the paper where its complexity is also discussed. 



Algorithm 1 Algorithm to find T min as stated in Lemma 3 



1: procedure [T] = Find_Tmin_One_Epoch (E, B lf B 2 , T upP er) 

2 . ymin ^_ Q) ymax ^_ T upper , T <- T max , E <- 

3: loop 

^. rp ^^min _|_ jimaxj 

5: E^T.g{B 1 /T,B 2 /T) 

6: if (E - E) > e then 

7: j-max ^_ rp ^ ess £ un j]; S f energy is used for transmission. In order to make E closer to E, we need to 
decrease T while transmitting exactly B\ and B2 bits. Reducing T max guarantees this operation in the next iterate.} 

8: else it(E-E)>e then 

9: jrmn ^_ ji |More than E units of energy is used for transmission. In order to make E closer to E, we need to 
increase T while transmitting exactly B\ and B2 bits. Raising T min guarantees this operation in the next iterate.} 

10: else 

11: return 

12: end if 

13: end loop 

14: end procedure 



The following extends the result of Lemma [3] to the case with two energy harvests. 
Lemma 4: Suppose Ei = Vz > 2 in the system model in Fig. [T] To minimize the overall transmission 
duration the sender finishes transmission to both users at the same time. 

Proof. The proof will make use of Lemmas [2] and [3] Take any rate pair (r n , r 2 i) for the I s * epoch resulting 
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in bn = rn£i, 6 2 i = r 2i£i bits being transmitted to the two users. This leaves &i 2 = B\ — bn, 622 = 
B 2 — 621 bits to be sent in the 2 nd epoch. By Lemma [3} total time to finish these remaining bits, which 
is T opt (2) = T min = min(max(Ti,T 2 )) will be minimized by setting T x = T 2 . More explicitly, 

^12 ^22 

T min = min(max(Ti, T 2 )) = min(max( — , — )) 
= mm(max( , )) 



5i-m£i 

ri2 



ri2 



^22 



rig ^ -Bi— rn^i 



»~22 B2— r21^1 

^2-^21^1 £12 > m|i 

r22 ron — riiA 



»*22 -B 2 -r 2 i^i 







A 























Infeasible 
region 



Bj-iiif;, 



Bi-fiiC, 
B;-r- : :, 



B,- r n';i 



decreases 



Fig. 4. Illustration of the proof of Lemma [4] 

T opt (2) is obtained by the rate pair satisfying r 12 /r 2 2 = (B 1 — r 11 ^ 1 )/(_B 2 — ?"2i£i) as point A on 
achievable rate region shown in Fig. |4j By Lemma [2| for the resulting two-epoch schedule to be optimal, 
a constant power and rate pair must have been used in the I s * epoch, and we have just proved that T\ = T 2 
for any constant choice in the I s * epoch. ■ 

Finally, we generalize the first result of Lemma |4[ to a general number of energy harvests. 

Lemma 5: Consider the system model with an arbitrary number of energy harvests described in Fig. 
[I] In a schedule that achieves T opt , the transmission to both users ends at the same time. 

Proof. The claim has been proved for k = 2 energy harvests, in Lemma |4j We will prove the general 
case by induction. Suppose that there are k energy harvests with the k th one at time tk, and the induction 
hypothesis holds, such that the optimal scheduler finishes transmission to both users at T opt (k). Now, 
consider adding a new energy harvest at time £fe+i- We have the following possible cases: 

1) tk+i > T opt (k): By the time the (k+l) st energy harvest arrives, the transmission has been completed, 
so by causality this energy harvest cannot help, and T opt (k + 1) = T opt (/c). 

2) t k+ i < T opt (k): In this case, the (k + l) st harvest will be used, to reduce the completion time. 
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Starting at time tk, the sequence of rate pairs will change from {(rij,r 2 j)} for i = to 
{(fii,f 2 i)} for i = l,...,k and (ruk+i), ^2(k+i)) f° r the newly added epoch. As in the proof of 
Lemma |4j whatever the number of bits allocated to the new epoch is, this rate pair will have a slope 
equal to the ratio of the number of bits remaining for this epoch. Hence, the bits will be terminated 
at some time T opt (A; + 1) < T opt (A;). ■ 

We are now ready to state the broadcast transmission scheduling problem as an optimization problem. 
From Section [TTJ for a given rate pair, there corresponds a unique power level P given by g(r 1; r 2 ) such 
that this rate pair is on the boundary of the rate region with power constraint P. The function g(ri,r 2 ) 
is strictly convex and continuously differentiable in r\ and r 2 . Using Lemmas [2] and [5} the problem can 
be written in terms of epoch rates. 

Given £> 1; B 2 , and the sequence {Ej}, supposing the problem is feasible (the total amount of energy is 
sufficient for transmitting the total number of bits), one can find an upperbound for the transmission 
completion time, T up in several ways. A simple one (which is possible when E\ is sufficient for 
transmitting the total number of bits), is to set the power so low such that only the first harvest is 
used to transmit all the bits. A much better upperbound will be obtained by the procedure that will be 
described within the initialization step of the FlowRight algorithm, in Section |IV} 

Given an upperbound for completion time, T up , we set k up equal to the index of the last energy harvest 
before this time, that is, k up = max{i : Yl)=i — ^ up }- An optimal solution will use at most k up harvests, 
and WLOG, remaining harvests can be ignored. Hence the problem reduces to finding T opt (/c up ): 

Problem 1: Transmission Time Minimization of Data Available at the Beginning on an Energy 
Harvesting Broadcast Channel: 

Minimize: T = T({(r li ,r 2 i)}i<i<k^p) 
subject to: r u ,r 2l > 0, 1 < i < k up 
< T < T up 

k k 

^2g(ru,r 2i )^<J2 E i (10) 

i 

for k = 1, 2, k* = ma,x{i : V] < T} 
fc* /,„ ^k' , ni / v^k* \+\ / k \ fe+1 



(11) 

i=l / i=l 



The set of constraints in (10) ensure that energy causality is respected. At any time during transmission, 
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the sender should have consumed at most the energy harvested up to that point, whereas by the end of 



transmission, it should have consumed all the harvested energies up to that instant. The constraint in ( fTTj ) 
ensures that all the bits of each user have been transmitted by the time T. Note that, by assigning nonzero 
values to all k up rates, one obtains a continuum of values of T that satisfy the constraints (note the 
( ) + used in the last constraint which sets the result to zero whenever the argument is negative), but the 
infimum of these, T opt , is the solution of the problem. We shall define n opt = min{i : J2)=i £j — ^ opt }> 
i.e., the index of the last harvest used by a solution that achieves T opt (It can easily be shown any such 
solution completely consumes all harvests from i to n opt .) 

This is not a standard convex optimization problem due to the objective appearing in the final equality 
constraint. Yet, we will establish that this minimization problem can be solved iteratively using an 
adaptation of the FlowRight algorithm [6|. Before moving on to the solution, we present our final 
observations in the optimal schedule in Theorems [I](2]l1 using the general result presented in Lemma [6] 

Theorem 1: In an optimal schedule, 

1) Powers assigned to epochs are monotonically nondecreasing, i.e., P\ < P<i < ... < P n ° P t. 

2) Energy consumed in any constant power band equals the total energy harvested within that band. 

3) The power assignment to epochs, P° pt = [Pi, P 2 , P n °pt], is unique. 
Proof. 

1) Suppose in an optimal solution we can find i s.t. Pi > Pi+\. From Lemma [TJ equalizing power 
over these epochs (this never violates causality as it corresponds to deferring the use of energy), 
one could find a rate pair with which more bits can be transmitted to each user. This contradicts 
the optimality of the original solution. 

2) Suppose that P { = P s ^ P s +i, s — m < i < s < n opt for some band of length m < s such 
that J2i=s-m E i Pt < Yfi= s -m E i- But b y Par t-1> P s +i > P s so a positive amount of energy (up 
to Ylri=s-m E i — J2i= s -m E i Pt an ^ s ) can be transferred from epoch s + 1 to s to equalize their 
powers, which by Lemma [T| allows rate pairs that send at least the same number of bits to each 
user. Keeping the rates and powers in the rest the same, the overall improved schedule is obtained, 
which contradicts the optimality of the original allocation. 

3) Suppose that there are two different optimal power allocation vectors, P A and P B , where P^ = P^ 
for k = l,2,..,i- 1 and P t A < P t B . From Part-1, power levels are monotonically nondecreasing 

2 Similar claims to those listed in Theorem [2] have been proven in 1111 through the observation that there is a cut-off level for the total 
power, below which the weaker user is assigned zero rate. 
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in the optimal schedule. In this case, if P^ for k > i stays constant, we have Y^k=i+i ^k^k < 
Ek=i+i P kZk> else 3 3 ■■ {P l A < Ptv 1 < 3 < n - i} and we have £tJ +1 P# fc < 
both contradicting Part- 2. ■ 
Lemma 6: Suppose the sender uses different rates for the stronger user in the intervals (ti, t*), (t*, t 2 ), 
such that T\ < r* < r 2 . Keeping powers levels and the number of bits transmitted to the stronger user in 
(ti, t 2 ) constant, a larger number of bits can be sent to the weaker user in (ti, r 2 ) by bringing the rates 
of the stronger user closer to each other if feasible. 

Proof. Consider the notation in Figure [2] and the case r n < r 12 . Keeping the avg. rate of user 1, 
rj., constant, set r n to r u , r 12 to r^ 2 s.t. r n < r u < r 12 < r 12 by transferring a certain amount of 
bits belonging to stronger user are transferred from the 2 nd slot to the I s *. This is feasible unless rn is 
already maximal for the given power level (i.e. r 21 = 0). Avg. rate of 2 nd user over the whole duration is 
increased from its original level, r 2 , to: 

f 2 = h 2 (P 1 ,r' 11 )f3 + h 2 (P 2 ,r' 12 )(l-^ 

> h 2 (Pi,r 11 )P + h 2 (P 2 ,r J3 )(l-P)=f 2 (12) 



( |12| ) follows from the fact that 

h 2 (P 1: r' n )P + h 2 (P 2 ,r' 12 )(l -P)- h 2 (P!, r n )P - h 2 (P 2 , r 12 )(l - P) > (13) 



for all /3 7^ {0, 1} (with equality achieved at P — 0, 1), unless r 21 = 0, as proved in App. IX-C In the 
remaining case, rn > r i2 , set r i2 < r^ 2 < r xl < rn, which is feasible unless unless r 22 = 0, and strictly 
improves the avg. rate for user 2. ■ 
Theorem 2: In an optimal schedule, 

1) the stronger user's rate is monotone nondecreasing, i.e., rn < ri 2 < ... < r ln o Pt ; 

2) if 7^ Tufor some < i < n opt , then r 2i = 0, i.e., if the stronger user's rate changes at the 
start of the (i + l) th epoch, the weaker user's rate was zero during the i th epoch; 

3) the weaker user's rate is monotone nondecreasing, i.e., r 2 \ < r 22 < ... < r 2n o P t; 

4) The vector of rate pairs, R opt = [(r°f ,r 2 f), (r°^ opt , r°^ opt )}, is unique. 
Proof. 

1) Suppose the rate of the stronger user decreases at some point, i.e., r^ > riu+n for some i. From 
Part-1 of Theorem [T] and Lemma |6j at least the same number of bits can be sent to each user (and 
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more to at least one) in epochs + 1) by assigning the strong user the average rate fi. 

2) Suppose that in an optimal schedule the weaker user's rate changes at the (i + l) th epoch and 
r 2i 7^ 0. If the rate of the stronger user changes, it can only increase, i.e., ru < rm + n, by Part-1. 
From Lemma |6} the schedule could only be improved by bringing ru and rm+i) closer to each 
other using the energy available for the weaker user at the i th epoch, if possible. Hence, the only 
reason why the rate of the stronger user can increase at the [i + l) th epoch is that there is no feasible 
energy available to equalize ru and ruj+i), which contradicts r 2i ^ 0. 

3) Suppose that in an optimal schedule r 2i > r 2 n + iy From Part-2, ru = ru i+ x\ = fi if r 2i ^ 0. By 
Part-1 of Theorem [I] and 2 nd property of the rate region, r 2i = h 2 (Pi,fi) < h 2 (P i+ i,fi) = r 2 (i+i) 
which contradicts initial rate assumption. 

4) To reach contradiction, suppose that there are two distinct optimal rate-pair vectors, R A and R B , 
where (V4,r 2 4 fc ) = (^ffc,^) for k = 1,2, — 1 and < r^. Using Part-3 of Theorem [T| we 
have r& = h 2 (Pi,r^) > h 2 (Pi,r^) > 0. From Part-2, rj- = r u < < rg V j e {i + 1, ...,n}. 
Hence, fewer bits will be transmitted by R A than R B , which contradicts the optimality of R A . ■ 

From Theorems [T] and [2j we conclude that the optimal schedule is unique (henceforth abbreviated as 
OPT.) The next section is devoted to the solution of Problem [T] 

IV. Optimal Offline Scheduling with the FlowRight Algorithm 

FlowRight stars from a feasible initial schedule, and progresses iteratively. Each iteration strictly 
improves the schedule (decreases T), which ultimately converges to the unique optimal T opt . 

Initialization: The energy consumed in each epoch is set precisely equal to the energy harvested at 
the beginning of that epoch. This schedule is feasible and P® = Ei/^. Given P®, j G {1,2,...}, one can 
assign rate pairs (r u , r 2i ) i G {1, 2, ...} on the achievable rate region boundary such that r u /r 2i = B\jB 2 . 
Let n up = argmin {r u = O,?^ = 0}. Algorithm |Ij presents a pseudo-code for this initialization. 

Algorithm 2 Initialization of FlowRight 

1: i <- 

2: while B 1 ^ 1 1 B 2 ^ do 
3: i++ 

4: Select (r° u , r° 2l ) such that: g(r° u ,r 2l ) EJ^ and r u /r° 2l ^ B 1 /B 2 
5: Bi Bi — T^Ci {Update remaining bits of I s * user.} 
6: B 2 B 2 — T 2i £,i {Update remaining bits of 2 nd user.} 

7: end while 

8: n up i {Set the initial number of epochs to be considered.} 
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After initialization, FlowRight performs local optimizations on pairs of epochs sequentially, i.e., on 
epochs (1,2), (2,3), (3,4), ... , until all epoch pairs are processed. This completes one iteration of the 
algorithm. Then, it continues with the next iteration, again performing local optimization on pairs of 
epochs at a time. The algorithm terminates after K iterations such that K = min {k : (T fc_1 — T k ) < e}, 
where T k is the transmission completion time at the end of k th iteration. 

Local optimization: Let E\ be the energy consumed at the i th epoch and n k < n up is the number of 
epochs used at the end of the k th iteration. Then, Ef = E{, % = 1,2, n np . Also, let be the number of 
bits transmitted to j th user at i th epoch at the end of k th iteration. Now, consider the epoch pair (i, i + 
i G {1, 2, ..,nfc — 1}. Local optimization aims to transmit the total b k = b^ 1 + ^7+i) Mts m me minimum 
amount of time by using the total available energy, i.e., E k ii+1 ^ = E\~ l + E\+l while respecting energy 
causality, i.e., using at most Y^ m =x ~ Em=i amount of energy in the i th epoch (see the term 
E^ ax . in Algorithm [3j). After a local optimization, we update the rate of the j th user in the i th epoch to 
r|j. These are the final values of the i th epoch rates as of the end of the k th iteration. We then reset the 
rate of the j th user in the (i + l) th epoch to r^Z+iy unless (i + 1) = n k . The (i + l) th epoch rates as of the 
k th iteration are finalized after the local optimization on the epoch pair (i + 1, i + 2) has been performed. 

Then, we continue with local optimization on the epoch pair (i + 1, % + 2). We proceed in this way to 
obtain (r^r^) for i = 1,2, ...,nfc. When the k th iteration is finished, we start from the beginning and 
update rates two epochs at a time similar to the above. It will be shown in Theorem [3] that transmission 
completion time strictly decreases after each iteration, and the number of epochs used, n k , is non-increasing 
from iteration to iteration. We terminate after K iterations, where K = min{/c : (T fc_1 — T k ) < e}. 
Algorithm [3] is the main pseudo-code for FlowRight. 

Details of the Local Optimization The local iteration step is described in Algorithm [4] which mainly 
checks whether it is possible to transmit b\ and b\ bits, (by respecting energy causality) in the minimum 
amount of time via equalizing power and rates in two epochs (it is optimal from Lemmas [T] and [6]), or 
not. If feasible, then min. transmission time is found by treating the epoch pair (i, i + as a single epoch 
with duration T uppcr = & + Then, Alg. [I] starting with the initial [T min = 0,T max = T uppcr ] domain 
will find the unique solution of local optimization. If power cannot be equalized, then we are at an energy 
causality boundary. In this case, Alg. [5} which can be implemented using the bisection method, is used 
to find the updated rate pairs (rj^rfj. The local optimization given in Alg. [5] exploits the structure of a 
locally optimized epoch pair (see Theorem [4] to work efficiently. The proof of this structure is omitted 
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Algorithm 3 FlowRight 



k <r- 0, T° <r- T U P 

repeat 

k++ 

for i = 1 to (nfc_i — 1) do 

°1 ^ °li + °l(i+l) 

<- min {(ELi ^ - . (Et 1 + 4+i))} 

= Fmd_Local_Optimal (E^ , E\~ x , E\~l , b\, b\) {Perform Local 

Optimization using Algorithm 13} 
end for 



if ==0&& 6 2 fc (nfc -i) ==0then 

nfe=nfc_i — 1 {Reduce the number of epochs for the next iteration.} 
else 

n k =n k -i 
end if 

Calculate_T (&T fc ) {Calculate current transmission completion time after one iteration over the epochs.} 
until (I*- 1 - T k ) < e 



for brevity but follows similar arguments on bisection method of Alg. [T] 



Algorithm 4 Algorithm to perform Local Optimization on the epoch pair + 

1: procedure [b k u , b k ^ iy b k t , b k 2 ^ iy E k , E^y = Find_Local_Optimal , E k ~ x , E^ 1 , b k , b\, ) 

3: T 1 min= Find_Tmin_One_Epoch {E^ ax . , b\, b\,^i + {Calculate T min to transmit &f and b\ bits using all the 

feasible energy for the i th epoch, i.e., E^ ax . at k th iteration.} 

4: if T min < & then 

5: {-^max- * s sufficient to transmit all the b x and 63 bits.} 

6: b k u <- , 6| <- 6f , 6^ x) <- , 6^ x) <- 0" {Update/Reset bits in the i t7 7(i + l) th epoch.} 

7: B* <r- g(r^, r^).T min , B^ 1 <- £$ i<+1) - {Update/Reset energy in the i th /(i + l) th epoch.} 

8: else 

9: T min = Find_Tmin_One_Epoch (B ( \ i+1) , 6f , + 

10: E k <r- g{b\£,i/T min , b k ^/T min ) {Update energy within the i th epoch} 

11: if E k < E k naXi then 

12: {We can equalize powers without violating energy causality} 

13: b k u <- b^/T mm , b\ t <- b^/T min {Update bits in the 1 st epoch.} 

14: b i(iii) <-bi- b k u , &2«!i) <~ & 2 - { Reset bits in the (* + i)"' epoch.} 

15: -E^ 1 <- - Bf {Reset energy in the (i + \) th epoch.} 

16: else 

17: {We can NOT equalize powers without violating energy causality. Maximum feasible energy should be allocated 
for the i th epoch.} 

18: E k <- E*^ , E^l <- B ( fe M+1) - {Update/Reset energy in the i th /(i + l) th epoch.} 

19: [r^,r^J=Find_Tmin_Two_Epoch(£ , f,B l fe 4 7 1 1 ,65 ; ,6|,^,^ + i) (see Algorithm p|. 

20: 6^ <- , 6^+1) <~ 6 i - 6 n ' 5 2* <~ r 2& - 6 2(m) <~ 5 2 - &2i {Update/Reset bits in the i tJ V(i + l) th epoch.} 

21: end if 

22: end if 

23: end procedure 



Next, we first establish that the cost function (that is, the completion time T) strictly decreases after 
each iteration of FlowRight, until it stops. We then show that the algorithm always stops. 
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Algorithm 5 Algorithm performing Local Optimization on the epochpair (i, i + i) when energy causality 
is met in the i th epoch. It is based on the observations in Theorem [4 



l 

2 
3 

4 
5 
6 
7 
8 
9 
10 

11 
12 
13 

14 
15 
16 
17 
18 
19 
20 
21 



22: 
23: 



24 
25 
26 
27 
28 



procedure [r llr r 2 i 
r^ in <- 0, 



Find_Tmin_Two_Epoch (E { , E i+1 ,b lr b 2 ,t; l , &+i 1 



'■g-^mm{/ ll (£ i /& ) 0) , iV&} 
rij <- r™ ax , r 1(i+1) <- 0, r 2l <- 0, r 2(l+1) <- 
{Check whether locally optimal rates of 1st user changes or not.} 
if /ii(£i/6,0) < then 

T<- Ox - r H ^)/r 1(i+1) 

r 2l «- h 2 {E l /£ !l ,r li ), r 2(l+1) «- h 2 (E i+1 /T,r 1{i+1) ) 
h = r 2i £i + r 2(i+1) T 
if (b 2 — b 2 ) > e then 

Find_Tmin_One_Epoch (-Ej+i , b\ — ru£i, b 2 — r2i£i r £i+i ) {Calculate minimum transmission time to transmit 
remaining bits within the (i + l) th epoch using Sj+i amount of energy} 
return 
end if 
end if 

{Locally optimal rates of 1st user are the same. Start Bisection Method.} 
loop 

ru = (r?T + r£ ax )/2 
r 1(i+1 ) <- r u 
T <- (h - rii&)/r 1(i+1 ) 
r 2i <- h 2 {E i /£i,r li ), r 2{i+1) 



h 2 {E i+ i/T, r 1{i+1} ) 



b 2 = r 2i £i + r 2 ( l+ i)T 
if (6 2 — b 2 ) > e then 



™ n guarantees 



r™ ax <r- ru {Less than b 2 bits are transmitted to the 2 nd user. In order to make those bits closer to b 2 , we need 
to decrease energy allocated to the 1 st user by increasing T while transmitting exactly bi bits. Reducing r™ ax 
guarantees this operation in the next iterate.} 
else if (b 2 — b 2 ) > e then 

r™ m ru {More than b 2 bits are transmitted to the 2 nd user. In order to make those bits closer to b 2 , we need to 
increase energy allocated to the I s * user by decreasing T while transmitting exactly 61 bits. Raising r 
this operation in the next iterate.} 
else 

return 
end if 
end loop 
end procedure 



Theorem 3: The following statements hold: 

1) As FlowRight runs, objective function T of P rob. ^strictly decreases after each iteration. Conversely, 
if T did not change after an iteration, then FlowRight stopped at the previous one. 

2) FlowRight stops, and returns a sequence {r^,r%°}. 

Proof. 

1) Suppose that we are on the k th iteration of the algorithm, and so far local optimizations have been 
performed on all epoch pairs up to (i — 1, i). During the (k — l) th iteration, local optimization on 
epochs (i + l,i + 2) has updated (r^T* , r 2(i+i)) wmi tne corresponding bits (bu^\, ^(i+i))- As 
the k th iteration progress, local optimization on the epochs (i — l,i) determines (rf^j,, r 2U-i)) 



17 

and resets (r^ , r^ -1 )- Suppose that this local optimization results in T mm before the end of the 
i th epoch. Hence in the rest of that epoch, rate pair changes to (0,0), i.e., a gap occurs. Let 
bj = bj° ^ + bj,7^1, j e {1,2} be the total bits to be transmitted in epochs + 1). From 
Lemma [2| by using a constant rate pair within the i th epoch, i.e., filling the gap within the epoch, at 
least the same number of bits can be transmitted to each user as in the original, slotted allocation. 
As the k th iteration progress, local optimization on the epochs + 1) determines (r^,^) and 
resets (^m+^j^fi+i)) by a t least trying to assign a constant rate pair in the i th epoch using the 
same amount of energy, hence reducing the number of bits to be transmitted in the (i + l) th epoch. 
Then, transmission surely ends before the end of the (i + l) th epoch, i.e., the gap moves from the 
i th epoch to the {% + l) th . Therefore, after local optimization, the total time to transmit b\ and b\ 
bits within the epoch pair + 1) reduces. Continuing this way, the gap propagates to the end of 
the k th iteration, hence the initial B\ and B 2 bits are transmitted by the time T k < T k ~ l . 
To prove the converse claim, suppose that after k th iteration, T did not change. Then, no gap has 
occurred during local optimizations; otherwise it would have propagated to the last epoch used and 
hence, reduced T. Therefore, performing further iterations can not create any gaps meaning that 
algorithm has indeed stopped at (k — l) st iteration. 
2) FlowRight initially starts from a feasible T > T opt , which is obviously lower bounded by T opt , 
the unique smallest completion time. From Part-1, T k = T({r^, r^}) is a strictly decreasing real 
sequence bounded below by T opt e 3ft, hence the iterations eventually stop. Therefore, the rate pairs 
{ r iii r 2i\ converge to some final value {r^°,r^}. ■ 

Now, we show in the following theorem that the schedule returned by FlowRight possesses the properties 
of OPT listed in Theorem Q] and |2l 
Theorem 4: When FlowRight stops, 

1) Powers are monotonically nondecreasing, i.e., P\ < P2 < ... < P n - 

2) Energy consumed during any constant power allocation band equals the total energy harvested in 
that band. 

3) The stronger user's rate is monotone nondecreasing, i.e., rn < r i2 < ... < r\ n , 

4) If the stronger user's rate changes at the {i + l) th epoch, the weaker user's rate is zero at the i th 
epoch, 

5) The weaker user's rate is monotone nondecreasing, i.e., r 2 i < r 22 < ■•• < r 2n . 



Proof. 

1) Suppose that we can find two epochs i and i + 1 with powers such that Pj > Pj+i. From Lemma [2j 
by bringing power levels closer to each other (this never violates energy causality), further local 
improvement on these epochs is ensured which contradicts the fact that FlowRight has stopped. 

2) Suppose that P, = P s ^ P s +i, s — m < i < s < n M for some band of length m < s such that 
V s Ef° < V s P°. Then, from Part-1, we have P s < P,+i. But we can transfer up to 
X^= s -m — J2i= s -m un i ts °f energy from epoch s + 1 to s, only improving the schedule (cf. 
Lemma [T]). This contradicts the assumption that FlowRight has stopped. 

3) Suppose 3 i G {l,..,?^ — 1} s.t. ru > ri( i+ i). From Part-1 and Lemma [6| by assigning the 
average rate f\ = (ru^i + ri(j + i)^ + i)/(^ + to the stronger user in epochs, more bits can be 
transmitted to the weaker user which means further local improvement on these epochs is ensured. 
This contradicts the assumption that FlowRight stopped. 

4) Suppose that the stronger user's rate changes at the (i + l) th epoch and r 2i ^ 0. The stronger 
user's rate can only increase, i.e., < ri(i + i), by Part-3. From Lemma |6j by bringing r u and 

closer to each other using the energy available for the weaker user at the i th epoch, overall 
transmission duration is decreased which contradicts the assumption that FlowRight stopped. 

5) Suppose that r 2 i > r^u+x)- From Part-5, ru = ri(j + i) since r 2 i ^ 0. From Part-3, the power is 
monotone increasing. Since h,2(P,r) is monotone increasing in P by the properties of rate region, 
r2(i+i) = h 2 {Pi+i,r li ) > h 2 {Pi,r u ) = r 2i , which contradicts r 2i > r 2 (i+i). ■ 

Theorem 5: The schedule returned by FlowRight is optimal, i.e., T({r^, r£°}) = T opt . 
Proof. Suppose that FlowRight stops and returns a schedule {r^,r^} = S' fr , with completion time 
nK°,r£})^T fr . 

As FlowRight respects feasibility, T fr can not be smaller than T opt . Suppose T fr > T opt . Consider the 
case that T opt is in the m th epoch and T fr is in the n th epoch with n > m. There must be a schedule 
{r°u,r° 2 f} = S opt that achieves T opt . Suppose that S opt and S ir are equal up to epoch s, which is the 
first time they differ either in terms of power level or rates, or both. Let us denote the power allocated in 
epoch s in S opt as P° pt and in S iT as P s fr . Consider the following. 

1) P s fr > P° pt : From Part-2 of Theorem [IJ all the harvested energies are consumed within any constant 
power band of 5 opt . Then, starting from epoch s when S opt consumes all the energy at the end of 
that constant power region, 5* fr would have consumed more energy than S opt by Part-1 of Theorem ffl 
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contradicting the fact that FlowRight always respects energy causality. 
2) P s fr < P° pt : Suppose that P s fr increases to P s f + m at some further epoch s + m before T opt . From 

Lemma [TJ by bringing power levels P^_ m _ 1 = Pf and P s f + m of the epoch pair (s + m — 1, s + m) 

closer to each other (this never violates energy causality), further local improvement on these epochs 

is ensured which contradicts the fact that FlowRight has stopped. 
Hence, S fr cannot have higher power level than S opt until T opt . Moreover, if power level of S ir becomes 
lower than that of S opt , then it should stay constant until T opt . These results are shown in the general case 
in Figure [5] Now, suppose that the first change occurs when P s fr = P° pt and the rate pairs, {r°£\ r^*} 



Same power schedule 
in this region 



Optimal 
FlowRight 



Same rate schedule Different rate schedule Different power & rate schedule 
in this region at k constant power bands at ' constant power bands 



Power 




opt 



Fig. 5. Illustration of the general case (in the proof of Theorem ||} that S b and S opt differ in power at k constant power bands and differs 
in both power and rate at I constant power bands. 



and {rf^rfj}, differ from each other in the general case. Consider the following. 



1) r°f < rf s : Since r°f = h 2 (P° p \ r°f ) > h 2 (P*,rf s ) > 0, rate of stronger user in the S opt should 
stay constant after epoch s by Theorem [2| Since T fr > T opt and rf s < rf u Vu G {s + 1, ...,n}, 
S ir should have transmitted more bits to stronger user than S opt , which contradicts the fact that 
FlowRight always respects bit feasibility, i.e., S ir transmits exactly the same number of bits to each 

user as S opt by the time T fr . 

2) r°f > rf s : We have rg = h 2 (P^ r , rfj > /i 2 (P° pt , r°f ) > 0; therefore, rf u = rf s \/u G {s+1, m} 
by Theorem 4j Then, bit feasibility requires rf ^Ei=a +1 t^j < Yl\=[^i r °i{iy wnere r °(£) > r iT 
is the rate of stronger user for the i th constant power band whereas ti is the duration of that 
band. Rearranging the terms we have rf < $^i=i7» r im> where ji = U/Y^^U £ (0, 1), 
Vi G {1, 2, ..,k + /}. Moreover, from Part-2 of TheoremjJ we have P fr **) > EtiS+i ^ P i> 
where p is the power of the constant power band (see Figjs]). Rearranging the terms we 

VES+i^ e (0,1), G {1,2,..,* + Z} and 



have P fr > /3ES + iK//?) p i where 




a ; 



E?=fc+i*i/ESfc+i*« e (0,1)- Now, let 6f and be the number of bits transmitted to 
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the 2 nd user from epoch s till T opt + t k+i+1 by S* and till T opt by S 10 ^, respectively. Then, we have 

k /k+l+1 \ /c+i 



i=l 

/fc+i+1 

E^ 

Vj=fc+1 
/fc+i+1 

\i=fc+l 
/fc+i+1 

E^ 

\i=fc+l 
/fc+i+1 

E^ 

\i=fc+l 

> £ * 

> £ * 

\i=k+i 

/k+l + 1 

= 



> 



> 



opt ' 



\i=fc+l / i=l 
fe fc+2 

E ai/i 2 (Pi, if ) + h 2 (P b , rf ) - E aM{Pi,r° 1{i) - , 

i=l i=l / 

k ( \ k+l \ 

E«« MJVft-MPi.r^j) +M^.»f)- E 



i=k+l 



J 



k + l 



h 2 (/3 e (oi/^.jf)- e «^2(p,^*) 



i=fe+l 

fe+; 



E («i/OT.»f)- E «^2(P,r^' 



i=k+l 



i=k+l 



k+l k+l \ 
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E 



\=k+l 



h^r^-h^r^) >0 



(14) 
(15) 
(16) 
(17) 

(18) 



>o 



Eq ( [14] ), Eq ( [15] ) and Eq ( [18] ) follows from the 2 property of the rate region while Eq ( [16] ) and 
Eq (fT7]) follows from the 3 rd . Hence, S ir transmits more bits to the 2 nd user than S opt , contradicting 



the fact that FlowRight always respects bit feasibility. 
Then, power allocation and rate pairs of S opt and S fr cannot differ, so S ir = S opt and T fr = T opt . 



V. Algorithm Complexity 

The core computational step in the algorithm is the local optimization (Algorithm [4]). This entails 
the solution of a nonlinear equation, which in our numerical computations has been done iteratively 
using the bisection method (detailed in Algorithms [T] and [5]). The exact number of iterations, hence the 
convergence rate depends on the selected tolerance level e, but in our experiments convergence time of a 
local computation is typically on the order of milliseconds. 

Maximum number of iterations of Algorithm [T] is proportional to log 2 (T upper /ei) whereas that of 
Algorithm [5] is proportional to log 2 {r^ / 62). Then, the worst case computation time of local optimization 
be Ciocai oc max{log 2 (T uppev /ei),log 2 (r^f /e 2 )}. For n epochs, n — 1 local optimizations are performed 
at each iteration. Hence, the worst case computation time of an iteration is C iter oc Ci oca i x n opt . 
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While FlowRight theoretically terminates when T does not change from one iteration to the next, in 
terms of implementation it would make sense to stop the iterations when the change is within some e of the 
average epoch size. The choice of e will determine the number of iterations and hence the linear scaling 
coefficient of complexity. In our extensive simulations based on randomized energy harvest amounts, epoch 
durations and bits to be transmitted, we observed that the number of iterations for achieving convergence 
sufficient for all practical purposes is on the order of (n opt ) 2 . Hence, the worst case computation time of 



FlowRight is oc Ci oca i x (n opt ) 3 which has polynomial time complexity in n 



opt 



VI. A Numerical Example 

Consider a two-user AWGN broadcast channel with bandwidth lOOKHz and noise power spectral density 
No = 10~ 13 Watts/Hz. Suppose that the path loss from sender to the stronger and weaker users are 70dB 
and 75dB, respectively. The sender needs to transmit 800Mbits to the stronger user and lOOMbits to the 
weaker user. Energy harvests of amounts [10, 10, 20, 40, 60, 70, 90, 180, 190, 100, 50, 30, 10] Joules 
arrive at [0, 2, 5, 7, 9, 10, 11, 13, 14, 15, 18, 20, 23] hours. 

FlowRight computes the final schedule shown in Fig. |6j The stopping criterion used in this example 
was successive iterations being within - Y^=i & where e = 10~ 9 . Running on a PC with Intel Core 2 
Duo CPU (2.26GHz) and 2GB RAM, the algorithm stopped after 59 iterations in 1.29 second. In the 
resulting optimal schedule, the last two energy harvests are not used. Note that in the final schedule, 
transmit powers remain constant during epochs, and are non-decreasing in time. 
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Fig. 6. A numerical example for the execution of FlowRight. The top figure represents the initial transmission completion time, T up =20.08 
hours, after the initialization phase and the final transmission completion time, T m in= 19.20 hours after the termination of the algorithm to 
transmit £>i=800Mbits and Bi= lOOMbits for the given energy harvest instants with the corresponding energy amounts. In the bottom figure, 
transmit powers are shown to be [1.11, 2.78, 5.56, 15.28, 23.30] mW for the durations [5, 2, 2, 4, 6.20] hours in the final schedule. The final 
schedule is {(rf?, r|?)}=[(1.6, 0),(L6, 0),(4.0, 0),(7.8, 0),(18.7, 0.6),(18.7, 0.6),(18.7, 0.6),(18.7, 4.1),(18.7, 4.1),(18.7, 4.1),(18.7, 4.1)]Kbps 
with durations as shown in the last figure. 
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VII. Conclusions and Future Work 

In this paper, we formulated and solved the offline transmission completion time minimization problem 
on an energy harvesting broadcast link. We have observed that, in the optimal solution, energy harvests 
may not necessarily be depleted at the end of each epoch, and could be deferred for later use. The schedule 
tries to "hurry up and be lazy" at the same time. The sender picks rates from the broadcast capacity region 
judiciously, such that it completes transmission to both users at the same time, T. In the optimal schedule, 
the powers are non-decreasing in time, so that transmission rate is highest toward the end. We have shown 
that the problem can be solved efficiently with a modification of the FlowRight algorithm. Our proposed 
algorithm starts with an upperbound on T and strictly improves it after every iteration or "pass" through 
the schedule, and stops when T converges to the optimal value. 

There are a number of directions for further work related to the problem presented in this paper. One 
of these is solving the offline minimization problem when data arrive during transmission, rather than 
being available in the beginning. Our preliminary work on this modification of the problem indicates that 
its solution has similar structural properties to the first problem, such as the powers being nondecreasing 
in time, and rates not changing between data arrival or energy harvest instants. Here, the optimal solution 
has more reason to be "lazy" in terms of transmission rate, as data will continue to come and it may be 
wise to save energy for future data arrivals. We believe that a further modified version of the iterative 
algorithm described in this paper solves this version of the problem. 

A second direction for further work is addressing the multiple-access version of this problem. There, 
energy harvests will be occurring at the senders, possibly at different points in time. Finding a distributed 
solution for that case may be a difficult yet interesting problem. 

Finally, another issue of interest is time-varying channel gain. The case of time-varying channel gain 
is interesting, and perhaps more meaningful to be setup as an online problem, rather than an offline 
problem, as channel gain variation is often difficult to predict (whereas energy harvesting times or packet 
formation times may be known ahead of time in some applications.) As the offline problem formulation 
has facilitated the analysis of the problem, going to an online formulation is arguably the most important 
challenge. While there are different ways to formulate the online problem, for example, as a dynamic 
control problem, our intuition is that approximate methods that leverage the offline formulation may be 
more tractable and insightful. 
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IX. Appendix 



A. Proof of Proposition [7J 

I s * and 2 nd order partial derivatives of hx(P,r) and h 2 (P,r) for the AWGN BC are as follows: 
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From ([2]) and ([3]), hi(P,r) and h 2 (P,r) are nonnegative. Monotonicity follows from ( fT9| ) and pO] ) as 
the signs of the first order partial derivatives of h^P, r) and h 2 (P,r) in their respective domains are 



always fixed. From pi) and p2] ), /ii(P, r) and /^(-P, r) are concave in power and, respectively rate, when 



the other parameter is held constant. The last property follows from (23) and rt24j). 



B. Proof of Lemma [7J 

Using ([8]), (|9]) can be written as 

fi{P) = h 2 {P± + (1 - /3)AP,ri)/3 + h 2 (P 2 - 0AP,ri)(l - /?) - h2(Pi,ru)0 - h 2 (P 2 , r 12 )(l - /3). 

The 2 nd order derivative of /i with respect to /3 is the following [^] 
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3 /i2 x and ft2„ represent the first order partial derivatives of /12 with respect to P and r, respectively. Second order partial derivatives of 
are represented by hi X!e , h,2 xy , ft.2„ x and h,2 yy 
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+0{h 2 „(P 1 + (1 - f3)AP,ri){-AP) 2 + h 2xy (Px + (1 - /3)AP,ri)(-AP)(r u - r 12 )} 
* . ' - — y ' 

<0 =0 

+f3{h 2y APi + (1 - P)AP, ri)(-AP)(m - r 12 ) + /i 2 „ v (Pi + (1 - 0)AP, ri)(r u - n 2 ) 2 } 

S v ' S v ' 

=0 <o 
+ (1 - p){h 2xx (P 2 - [3AP,fi)(-AP) 2 + h 2xy (P 2 - )9AP,fi)(-AP)(rn - n 2 )} 

<o =0 
+ (1 - P){h 2yx {P 2 - ^AP,ri)(-AP)(rn - r 12 ) + h 2yB (P 2 - ri)(m - r 12 ) 2 } < (26) 

N ^ / S y / 

=0 <0 



According to the properties of the rate region (l)-(4), (26) always holds. Hence f\ is concave in (3. 



C. Proof of Lemma [6] 

Substituting r n = r u + (1 — /3)Ar and r' 12 = r i2 — /3Ar in to ( fT3| ), we have the following. 



HP) = h 2 (Pi,r n + (1 - P)Ar)/3 + h 2 (P 2 , r n - /3Ar)(l - p) - h 2 {Pi,r u )p - h 2 (P 2 ,r 12 )(l - (3). 
The 2 nd order derivative of / 2 with respect to (3 is the following 

= 2{h 2y (P 1 ,m + (1 - (3)Ar){-Ar) - h 2y {P 2 ,r 12 - pAr){-Ar)} 
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According to the properties of the rate region (l)-(4), (27) always holds. Hence / 2 is concave in (3. 
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